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One-dimensional scattering problem admitting a complex, VT— symmetric short- 
range potential V(rr) is considered. Using a Runge-Kutta-discretized version of 
Schrodinger equation we derive the formulae for the reflection and transmission coef- 
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1 Introduction 



Standard textbooks describe the stationary one-dimensional motion of a quantum 
particle in a real potential well V(x) by the ordinary differential Schrodinger equation 

ip(x) — Eip(x) , x G (—00, 00) (f) 

which may be considered and solved in the bound-state regime at E < V(oo) < +00 
or in the scattering regime with, say, E = k 2 > V(oo) = 0. In this way one either 
employs the boundary conditions ^(±00) = and determines the spectrum of bound 
states or, alternatively, switches to the different boundary conditions, say, 

f Ae iKX + Be~ iKX , x < -1 , 
ip(x) = < (2) 

[ Ce lKX , 

Under the conventional choice of A = 1 the latter problem specifies the reflection 
and transmission coefficients B and C, respectively [I]. 

The conventional approach to the quantum bound state problem has recently 
been, fairly unexpectedly, generalized to many unconventional and manifestly non- 
Hermitian Hamiltonians H 7^ which are merely quasi-Hermitian, i.e., which are 
Hermitian only in the sense of an identity = 6 H 6™ 1 which contains a nontrivial 
"metric" operator = > as introduced, e.g., in ref. [2]. The key ideas and 
sources of the latter new development in Quantum Mechanics incorporate the so 
called VI — symmetry of the Hamiltonians and have been summarized in the very 
fresh review by Carl Bender [3]. This text may be complemented by a sample [4] of 
the dedicated conference proceedings. 

In this context we intend to pay attention to a very simple VT— symmetric scat- 
tering model where 

V(x) = Z{x) + i Y(x) , Z{—x) = Z(x) = real , Y(—x) = —Y(x) = real 



dx 2 



+ V(x) 
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and where the ordinary differential equation (f) is replaced by its Runge-Kutta- 
discretized, difference-equation representation 

- 2«s t ) + ^) + v{xtmxk) = (3) 

with 

Xk = kh, h>0, k = 0,±l,... 
as employed, in the context of the bound-state problem, in refs. [5]. 

2 Runge-Kutta scattering 

Once we assume, for the sake of simplicity, that the potential in eq. (3) vanishes 
beyond certain distance from the origin, 

V(x ±3 ) = j M.M • 1 

we may abbreviate ipk = ip(xk), Vk — h 2 V(xt) and 2 cosy? = 2 — h 2 E in eq. (3), 

-ipk-i + (2 cos ip + V k ) i>k ~ i>k+i = . (4) 

In the region of \k\ > M with vanishing potential 14 = the two independent 
solutions of our difference Schrodinger eq. (4) are easily found, via a suitable ansatz, 
as elementary functions of the new "energy" variable ip, 

ipk = const ■ g k =>- q — q± — exp(±i (p) . 

This enables us to replace the standard boundary conditions (2) by their discrete 
scattering version 

( Ae im^ + Be -im<p ^ m < _( M _ ^ ^ 

tp(x m ) = < (5) 

[ Ce tmtp , m > M - 1 

with a conventional choice of A = 1. 
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Two comments may be added here. Firstly, one notices that the condition of the 
reality of the new energy variable tp imposes the constraint upon the original energy 
itself, -2 < 2 - h 2 E < 2, i.e., E e (0,4//i 2 ). At any finite choice of the lattice 
step h > this inequality is intuitively reminiscent of the spectra in relativistic 
quantum systems. Via an explicit display of the higher 0(h 4 ) corrections in eq. (3), 
this connection has been given a more quantitative interpretation in ref . [6] . 

The second eligible way of dealing with the uncertainty represented by the 0{h A ) 
discrepancy between the difference- and differential-operator representation of the 
Schrodinger's kinetic energy is more standard and lies in its disappearance in the 
limit h — > 0. This is a purely numerical recipe known as the Runge-Kutta method 
[7] . In the present context of scattering one has to keep in mind that the two "small" 
parameters h and 1/M may and, in order to achieve the quickest convergence, should 
be chosen and varied independently. 

3 The matching method of solution 

3.1 The simplest model of the scattering with M = 1 

Once we are given the boundary conditions (5) the process of the construction of 
the solutions is straightforward. Let us first illustrate its key technical ingredients 
on the model with the first nontrivial choice of the cutoff M — 1. In this case our 
difference Schrodinger eq. (4) degenerates to the mere three nontrivial relations, 

-■0-2 + 2 COS if - IpQ^ = 

-V>-i + (2 cos if + Z ) - 0i = (6) 
+ 2 cos up 01 -02 = 

where we may insert, from eq. (5), 

f^e^ + Be'^, 4" ) = l + #, 4 +) =C, ^ = Ce l * (7) 
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and where we have to demand, subsequently, 

4-) = i + J B = 4+) = c = ^ , 

-e^t - Be** + (2 cosy? + Z Q ) C - C e i(p = . 
Thus, at an arbitrary "energy" (p one identifies B — C — 1 and gets the solution 

= 2i sin y? ^ Z 



2i sin </? — Z 2i sin </? — Z 

Of course, as long as we deal just with the real "interaction term" Zq, our M = 1 
toy problem remains Hermitian since no VT— symmetry has entered the scene yet. 

3.2 VT— symmetry and the scattering at M = 2 

In the next, M = 2 version of our model we have to insert the four known quantities 

i{>_ 2 = e- 2iv + B e 2iv , V-i = e~ i(p + B e iifi , fa = C e iv , fa = Ce 2itfi 
in the triplet of relations 

-^_ 2 + (2 cos y> + Z_i - i y_i) V-i - 4 _) = 

+ (2 cos + Z ) - fa = (9) 
+ (2 cos + Z_i + i - ^ 2 = o 

where the three symbols ipo, ipQ~^ and ipQ + ^ defined by these respective equations 
should represent the same quantity and must be equal to each other, therefore. 
Having this in mind we introduce C,^ = 1 + B and £^ = C and decompose 



This enables us eliminate 
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and eq. (9) becomes reduced to the pair of conditions, 

l + B + V-^-! = C + V 1 ^ 1 =Vo, (1Q) 
-■0-1 + (2 cos v? + Z ) i/j - 0i = 

They lead to the two-dimensional linear algebraic problem which defines the reflection 
and transmission coefficients B and C at any input energy ip. The same conclusion 
applies to all the models with the larger M. 



4 The matrix- inversion method of solution 



Let us now re-write our difference Schrodinger eq. (4) as a doubly infinite system of 
linear algebraic equations 



5_i -1 
-1 S -1 
-1 Sx 



V 



\ / : 

00 

) \ '■ 



\ 



= 0, 



/ 



where 



Sk 



s*- k ) = 



(12) 



2 cos (p + Z k + i Y" fc sign k , |fc|<M, 
2 cosy^, | A; | > M 

and where the majority of the elements of the "eigenvector" are prescribed, in ad- 
vance, by the boundary conditions (5). Once we denote all of them by a different 
symbol, 



Ae imv + Be -imv = ^H > 



m < — (M — 1), 
m > M — 1 , 



(13) 
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we may reduce eq. (11) to a finite-dimensional and tridiagonal non-square-matrix 
problem 



-1 S 



(M-l) 



-l sz -l 

-1 S -1 

-1 S! -1 



\ / *(-) 

S-M 

S-(M-l) 
■0-(M-2) 



^M-2 
SM-1 



V 



-1 5, 



(M-l) 



v * 



(+) 

M 



(14) 



or, better, to a non-homogeneous system of 2M — 1 equations 



/ c(-) \ ( A-) \ 

S-M 



£-(M-l) 
■0-(iW-2) 

\ SM-1 / 











(15) 



where the (2M — 1)— dimensional square-matrix of the system can be partitioned as 
follows, 



T = 




(16) 



(M-l) J 



Whenever this matrix proves non-singular, it may assigned the inverse matrix R=T 1 , 
the knowledge of which enables us to re- write eq. (15), in the same partitioning, as 



6 



follows, 



( A-) \ 

S-(M-l) 

A+) 

\ SM-l / 



( A-) \ 

S-M 



R 



o 



W-{M-2) 



V 



(17) 



/ 



■0M-2 

In the next step we deduce that the matrix R has the following partitioned form 

\ 





( * 

a* 


P 





R = 


u 


Q 






{ P 







We may summarize that in the light of the overall partitioned structure of eq. (17), 
the knowledge of the (2M — 3)— dimensional submatrix Q as well as of the two 
(2M — 3)— dimensional row vectors i? and w T (where T denotes transposition) is 
entirely redundant. Moreover, the knowledge of the other two column vectors u and 
v only helps us to eliminate the "wavef unction" components ■ifj_( M _ 2 ), ip-(M-3), • • • ? 
■0M-3, i>M-2- I n this sense, equation (15) degenerates to the mere two scalar relations 



A-) -rA+) 

S-(M-l) " S-(M) 



: 0. 



(18) 



Mt; _ _ 
SM-l h"r,-M u SM 

Once we insert the explicit definitions from eq. (13) we get the final pair of linear 
equations 

e -HM-l)<p + Be i(M-l) v _ a * ^ e -iM V + Be iM^ -C (3e iM * = 0, 

which are solved by the elimination of 

C 



(19) 



B = -e 



-2iMtp 



+ 



and, subsequently, of 



C 



{eT* - a) 



smip 



This is our present main result. 



(20) 



(21) 
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5 Coefficients a and f3 

Our final scattering-determining formulae (20) and (21) indicate that the complex 
coefficient a and the real coefficient (3 carry all the "dynamical input" information. 
At any given energy parameter <p these matrix elements are, by construction, rational 
functions of our 2M — 1 real coupling constants Z , Z±, . . . , Zm-i and Y±, . . . , Ym-i- 
In particular, j3 is equal to 1/ det T and a has the same denominator of course. An 
explicit algebraic determination of the determinant det T and of the numerator (say, 
7) of a is less easy. Let us illustrate this assertion on a few examples. 

5.1 M = 2 once more 

det T = Z Zj 2 — 2 Zj + Yi 2 Z 

Re7 = Z Z 1 -\ 
Im7 = -Z Y t 

5.2 M = 3 

det T — Zg Z \ ^ Z2 ^ — 2 Z\ Z2 — 2 Z2 ^ ~h 

+ y_/ Zq Z2 + 2 + y# Zq Zi — 2 y# -\- Y2 Yi Zq + 2 y^ y# + 

Re7 = Z - 2 Z g + Y, 2 ^ ^ - Z t Z + 1 

Im 7 = -Z Z 1 2 Y 2 + 2Z 1 Y 2 - Y t 2 Z Y 2 - Y 2 Z 

5.3 M = 4 

The growth of complexity of the formulae occurs already at the dimension as low 
as M = 4. The determinant det T and the real and imaginary parts of 7 are them 
represented by the sums of 15 and 14 and 32 products of couplings, respectively. 
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A simplification is only encountered in the weak coupling regime where one finds 
just two terms in the determinant which are linear in the couplings, 

det T = -2 Z 3 - 2 Z t + . . . 

being followed by the 10 triple-product terms, 

• • • + Y i 2 Zo + 4 Zi Z2 Z 3 — 4 Zi Y2 Y 3 + Z 2 Zq + 2 Y 3 2 Z 2 + 2 Yj Zq Y 3 -\- 

+2 ij/ + Z 3 2 Zq + Y 3 2 Zq + 2 Z 2 Z2 + . . . 

etc. Similarly, we may decompose, in the even-number products, 

Re 7 = -1 + 2 Z 2 Z 3 + Z Z t + Z Z 3 + 2 Z 2 Z t + . . . 

and continue 

• • • H — 2 Zo Zi Z2 Z$ + 2 Zo Yi Y2 Z 3 — 
— Z2 Y 2 Zo — 2 -Zg 2 Z\ Z 3 — Z2 Z 2 Zq — 2 Yg 2 Z\ Z 3 + . . . 

etc, plus 

Im 7 = — 2 Z 2 Y 3 + 2 Y 2 Zt — Z Y t — Z Y 3 — . . . 
with a continuation 

... — Fg Fj 2 Z# — Yg Z_/ 2 ^+ 

+2 Yg ^ Y5 + 2 Zg Z 2 Zi Y 3 — 2Z Y 1 Y 2 Y 3 + . . . 

etc. Symbolic manipulations on a computer should be employed at all the higher 
dimensions M > 4 in general. 

6 Discussion 

The main inspiration of the activities and attention paid to the VT— symmetry 
originates from the pioneering 1998 letter by Bender and Boettcher [8] where the 
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operator V meant parity and where the (antilinear) T represented time reversal. Its 
authors argued that the complex model V(x) = x 2 (ix) 5 seems to possess the purely 
real bound-state spectrum at all the exponents 5 > 0. After a rigorous mathematical 
proof of this conjecture by Dorey, Duncan, Tateo and Shin [9] and after the (crucial) 
clarification of the existence of a nontrivial, "physical" Hilbert space H where the 
Hamiltonian remains self- adjoint [2, 10, 11, 12], the bound-state version of eq. (1) 
may be considered more or less well understood, especially after it has been clarified 
that the physics-inspired concept of VT — symmetry of a Hamiltonian H should in 
fact be understood, in the language of mathematics, as a V— pseudo-Hermiticity of 
H specified by the property # f = VHV~ 1 [10, 11, 13]. 

In the spirit of the latter generalization, current literature abounds in the studies 
of the potentials which are analytically continued [14], singular and multisheeted 
[15], multidimensional [16], manybody [17], relativistic [18], supersymmetric [19] and 
channel-coupling [20] . Among all these developments, a comparatively small number 
of papers has been devoted to the problem of the scattering. For a sample one might 
recollect the key reviews [21] and various Kleefeld's conceptual conjectures [22] as well 
as a very explicit study of the scattering by the separable VT — symmetric potentials 
of rank one [23] or by the rectangular or reflectionless barriers [24], or the motion 
considered along the so called tobogganic (i.e., complex and topological^/ nontrivial) 
integration contours [25]. In this context our present difference-equation-based study 
may be understood just as another attempt to fill the gap. 

Technically we felt inspired by our old Runge-Kutta-type discretization of the 
VT— symmetric Schrodinger equations [5] as well as by our recent chain-model ap- 
proximations of bound states in a finite-dimensional Hilbert space [26]. In a certain 
unification of these two approaches we succeeded here in showing that there exists 
a close formal parallelism between the description of the (one-dimensional, Runge- 
Kutta-approximated) scattering by a real (i.e., Hermitian) potentials and by their 
complex, VT— symmetric generalizations. We showed that in both these contexts, 
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the definition of the transmission and reflection coefficients has the same form [cf., 
e.g., eq. (21)], with all the differences represented by the differences in the form of 
the "dynamical input" information. It has been shown to be encoded, in both the 
Hermitian and non-Hermitian cases, in the two functions a and (5 of the lattice po- 
tentials, with the vanishing or non- vanishing coefficients Y k , respectively (cf. a few 
samples of the concrete form of a and f3 in section 5). 

On the level of physics we would like to emphasize that one of the main dis- 
tinguishing features of the scattering problem in VT— symmetric quantum mechan- 
ics lies in the manifest asymmetry between the "in" and "out" states [22]. In its 
present solvable exemplification we showed that such an asymmetry is merely formal 
and that the problem remains tractable by the standard, non-matching and non- 
recurrent techniques of linear algebra. A key to the success proved to lie in the 
partitioning of the Schrodinger equation which enabled us to separate its essential 
and inessential components and to reduce the construction of the amplitudes to the 
mere two-dimensional matrix inversion [cf. eq. (18)] where all the dynamical input 
is represented by the four corners of the inverse matrix R = T" 1 [cf. the definition 
(16)]. 

We believe that the merits of the present discrete model were not exhausted by 
its present short analysis and that its further study might throw new light, e.g., on 
the non-Hermitian versions of the inverse problem of scattering. 
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